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Abstract. The double exchange (DE) model with quantum local spins S is studied; 
an equation of motion approach is used and decoupling approximations analogous to 
Hubbard's are made. Our approximate one-electron Green function G is exact in the 
atomic limit of zero bandwidth for all S and band filling n, and as n — > reduces 
to a dynamical coherent potential approximation (CPA) due to Kubo; we regard 
our approximation as a many-body generalisation of Kubo's CPA. G is calculated 
self-consistently for general S in the paramagnetic state and for S = 1/2 in a state 
of arbitrary magnetization. The electronic structure is investigated and four bands 
per spin are obtained centred on the atomic limit peaks of the spectral function. A 
resistivity formula appropriate to the model is derived from the Kubo formula and the 
paramagnetic state resistivity p is calculated; insulating states are correctly obtained 
at n — and n — 1 for strong Hund coupling. Our prediction for p is much too small 
to be consistent with experiments on manganites so we agree with Millis et al that the 
bare DE model is inadequate. We show that the agreement with experiment obtained 
by Furukawa is due to his use of an unphysical density of states. 



PACS numbers: 75.20.Hr, 72.90. +y, 71.28.+d 



1. Introduction 



Manganite compounds exhibiting colossal magnetoresistance (CMR) are of the form 
Lax^^D^MnOa with D divalent, e.g. Ca, Sr, Ba. As the doping x and temperature T are 
varied a rich variety of phases are observed, as discussed by Ramirez [p]] . Recently there 
has been a lot of interest in these compounds with x > 0.15 owing to their interesting 
magnetotransport properties in this regime: as T is decreased they undergo a transition 
to ferromagnetic order, and near the Curie temperature Tq the T-dependence of the 
resistivity p changes from insulating (dp/dT < for T > Tq) to metallic (dp/dT > 
for T < Tq), with a strong peak in p at the crossover. The application of a strong (~ 5T) 
magnetic field substantially reduces this peak in p and shifts it to higher temperature, 
giving rise to a very large negative magnetoresistance. The physical processes causing 
this behaviour have been the subject of much discussion. 
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The simplest model proposed for the CMR compounds — the one that we will study 
in this paper — is Zener's double exchange (DE) model with Hamiltonian 

H = X; UA^ia -J y ZSi'*i = H + H 1 . (1) 

ija i 

Here i and j refer to sites of the (approximately) simple cubic lattice of Mn atoms, 
c ja ( C L) is a c-spin conduction electron annihilation (creation) operator, Si is a local 
spin operator, &i is a conduction electron spin operator, is the hopping integral with 
discrete Fourier transform tr, and J > is the Hund's rule coupling constant. The 
number of conduction electrons per atom n is assumed to be given by n = 1 — x. 
Physically, the relevant electrons are those coming from the Mn atoms' 3d shells, which 
contain four electrons per site in the undoped compounds and are split by the cubic 
crystal field into triply degenerate t 2g levels and higher energy doubly degenerate e g 
levels. Strong Hund coupling attempts to align all electron spins on a site, so the t 2g 
electrons are treated as localized S — 3/2 spins while the conduction band is formed 
from the e g states. The main physical effects neglected by H are the double degeneracy 
of the e g conduction band, impurity atom (D) disorder scattering, and coupling to the 
lattice degrees of freedom. 

Furukawa || has studied the infinite dimensional limit of H for S = oo using 
dynamical mean field theory, and has concluded that the DE model's predictions for p 
in the paramagnetic state are compatible with experiment. Millis et al [f|, |5j however 
have claimed that p predicted by the DE model is much smaller than that measured, and 
that to get agreement with experiment dynamical Jahn- Teller phonon coupling must be 
included in H . Experiments show that phonon coupling is important, with for example 
a large shift in Tq observed upon replacing some of the O atoms with a different isotope 
||, so Furukawa's claim is puzzling. In this paper we will study the one-electron local 
Green function G and calculate p in an attempt to reconcile the results of Millis and 
Furukawa. This work has been briefly summarised elsewhere [7], [8| and is discussed in 
more detail by Green ||. 

Our starting point is Kubo's calculation based on a dynamical coherent potential 
approximation (CPA). For finite local spins S dynamical scattering processes may occur 
in which local spins and conduction electrons exchange angular momentum, whereas in 
the classical S — > oo limit taken by most authors the local spins are rigid and H is 
a one-electron Hamiltonian with spin dependent diagonal disorder. Kubo's CPA is an 
extension of the familiar alloy CPA which takes these dynamical processes into account 
in a local approximation. Since it is a one-electron theory Kubo's approximation is only 
valid in the low-density n — > limit, but in this limit the behaviour of the spectral 
function is qualitatively correct with bands with the correct weights forming about the 
two atomic limit (t^ — * 0) peaks as the hopping is switched on. If J ^> tij double 
occupation of a site is forbidden so that at half-filling (n = 1) the system should be a 
Mott insulator. This is not the case in Kubo's CPA where the Fermi level lies within the 
lower band. The correct behaviour will be obtained in an approximation which becomes 
exact in the atomic limit for all filling, so we are looking for a many-body extension of 
Kubo's CPA, valid for all n, which reduces to Kubo's CPA as n — > and to the correct 
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atomic limit as £y — > 0, for all n. 

It is difficult to extend the usual CPA method to the many-body case. Instead we 
return to the original approach of Hubbard [|1(]] in which he applied the equation of 
motion method to calculate the one-electron local Green function G for the Hubbard 
model. His decoupling approximation was motivated by the alloy analogy in which 
electrons of one spin are considered as frozen on atomic sites. Hubbard's 'scattering 
correction' is equivalent to a CPA treatment of the alloy analogy [|n]] and his 'resonance 
broadening correction' was an attempt to restore some dynamics to the frozen electrons. 
The idea of j spins being frozen in the calculation of the f spin Green function is 
introduced in the equation of motion method by neglecting commutators of the kinetic 
part of the Hamiltonian Hq with l-spin occupation numbers. 

In the derivation of our approximation to G for the DE model we will make 
approximations analogous to those used by Hubbard to obtain the scattering correction; 
terms corresponding to his resonance broadening correction will be neglected. Our 
method represents a considerable extension of Hubbard's owing to the more complicated 
form of the interaction term of the DE model, which for instance allows electrons to 
change spin via exchange of angular momentum with the local spins. This effect couples 
the equations for G^ and Gi which may perhaps be regarded as including some resonance 
broadening effects. Our choice of approximations will be guided by the requirement that 
we recover Kubo's CPA as n —* and the correct atomic limit as Uj — > 0. Owing to the 
spin symmetry of H we only need to derive an equation for G^, and the equation for Gj_ 
follows immediately. 

We calculate the atomic limit Green function g exactly in section |2| and derive our 
many-body CPA equation for the Green function G in section [5| The CPA spectral 
function is studied in section £|. In the limit of infinite J and n — > 1, Kubo [12] has 
introduced a CPA treatment of holes and we compare our results with this in section 
section |5|. A formula for the resistivity p of the zero field paramagnetic state is derived 
in section ||] and in section |7] p is calculated for various approximations to the density of 
states (DOS). A summary is given in section [8]. 



2. Atomic limit 



Since we require our approximation for G to be exact in the atomic limit we must 
first derive G in this limit. We define the retarded Green function for operators A 
and B with no explicit time- dependence by ((A;B)) t = —i6 + (t)({A(t),B}) and its 
Fourier transform by ((A;B)) e = /f^ dt e iet ((A ;B)) t . Here 9 + {x) = 1 for x > and 
otherwise, e is restricted to the upper half of the complex plane, and [, ] and {, } are the 
commutator and anticommutator respectively. The equation of motion for the latter 
Green function is 

e((A;B)) e = ({A,B}) + {([A,H];B)) e . (2) 

The cx-spin one-electron Green function is given by G l J(e) = ((c ia ; c] .)) e , and we define 
G G = G™ and G^ a = N" 1 J2r- ex P(i^ ' Rij)Ga where Rij is the vector from the j th to 
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the i th lattice site and N is the number of sites. The occupation number operators are 

±, such that ntr = n 4 „ and 



n 



11: 



11: 



c\ a c ia and we introduce the notation nf a , with a 



We also define the total occupation number n, = n»j + and 



1 



11: 



We denote the Green function G„ in the atomic limit 



by g 



In the atomic limit the eigenstates and energy eigenvalues of H are completely 
defined in terms of those of the single-site Hamiltonian. In order of increasing energy 
these consist of 2(5 + 1/2) + 1 = 25 + 2 singly occupied states with electron spin and 
local spin parallel having energy —JS/2, 25 + 1 unoccupied states and 25 + 1 doubly 
occupied states with energy 0, and 2(5 — 1/2) + 1 = 25 singly occupied states with 
electron spin and local spin antiparallel having energy J(S + l)/2. is easily obtained 
using the equation of motion method. We define = ((n a Cp cj)), = ((n a 5 2 Cp cj)), 
and = ((n a 5~C|; cj)), dropping site indices since the site referred to is always the 
same. S and a are the raising/ lowering operators for the local spins and conduction 
electrons respectively. We start with g a and write equations of motion for undetermined 
Green functions until the system of equations closes, obtaining 

J 



nr) — — 



(e)+t? 



S «(e) 



SZn l 



c- — Ja 

a b o ) 



e + 1\ 



JS{S + 1) 



-9" 



(3 a) 
(36) 



Now fl>f = J2 a =±9^, an d solving (|3^) and (^) for g a we obtain 
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25+1 



((5 + 5> ; - 5-a+) , ((5-5 z )(l-n i )-5-a+) 



+ 



e + J(5 + l)/2 
((5 + l + 5 2 )(l-n i )+5- 



+ 



J(5 + l)/2 
((5+1-5 2 



In i 



b-(T + 



(4) 

e + J5/2 e-J5/2 

It should be noted that for finite 5 there are four peaks in the spectral function at 
energies ±J(5 + l)/2 and ± J5/2; for 5 — > oo the upper and lower pairs of peaks merge 
leaving the familiar double-peaked spectral function. To obtain <?j,(e) from (|j) we make 
the changes n a \— > 5 Z i— > — 5 2 , 5 1 * 1 i— > 5 T , (T 1 * 1 h- > (T T . The weight summed over 
spin in the lowest band is [nb — 2(5 ■ a)) /{2b + 1), and since (5 ■ a) — > n5/2 as the 
temperature T — > for sufficiently large J the peak at — J(S + l)/2 will have very 
little weight at reasonable temperatures. Hence it will often be a good approximation 
to neglect the lowest band. The total weight summed over spin in the lowest two bands 
is (25 + 2 — n)/(2S + 1), so it is clear that a ty ^ theory that becomes exact as — > 
will give an insulator at half-filling. 



3. CPA Green function 



Our general strategy in this section follows that of Hubbard — we will retain occupation 
numbers as operators rather than replacing them with their averages, but their 
commutators with Hq will be neglected where appropriate, corresponding to the frozen 
electron approximation of the alloy analogy. Terms which are not multiplied by tij must 
be treated exactly in order to obtain the correct atomic limit. 



We split GV into two components, GV = J2 a =± where G l J a (e) = ((nfc^ ; cL)) e , 
corresponding to propagation through singly (a = — ) and doubly (a = +) occupied 
sites. Neglecting [nf^, H ] the equation of motion for each component is then given by 

eGf (e) » {nl) L + E^(e)) + J> ((K 6 *! ' c 1t)) 6 

\ k / k 

~(sr(e)+T^(e)) (5) 

where we have introduced the Green functions ^"(e) = ((nfS^c^ ; c]^)) e and T| ja (e) = 
{{nfS~ ; c]|)) e and the notation 5 A = A — (A) for any operator A. The second term 
on the right-hand side of (||) is Hubbard's scattering correction in which the deviation 
of nf^ from its average is accounted for. The last term of (|^) is more complicated than 
in the case of the Hubbard model, containing the as yet undetermined Green functions 
5*| and 7\. These correspond respectively to propagation of the electron as an j-spin 
and, following spin-flip scattering from a local spin, as a j-spin; the presence of T-f will 
couple the equations for G^ and Gj. 

We first treat the scattering correction, splitting the relevant Green function into 
two components, {{Snf^c k ^ ; cj|)) e = J2p=±((8nf^n k c k ^ > c ji))f ^ ^ s assumed that tu = 0, 
so from (|5|) it may be seen that this Green function is needed only for i ^ k. The 
equations of motion are, for i ^ k, 



<foi>ifcC fcT ; ct T )) « (Sn^Snfy 5 jk + J2hi ((Snf^c^ ; ct T }} 



J 
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((5nf L nt(S z k c k1 + S k c kl );c^)) e (6) 



where [nf|,i?o] and all commutators involving 8nf^ have been neglected. This is 
consistent with the strategy stated above. As a further approximation the first term 
on the right-hand side of @, a two-site correlation function, is dropped and we set 
rifci ~ ( n ki) i n the second term, which corresponds to neglecting a second scattering 
correction. 

The system of equations for the scattering correction is now closed apart from the 
last Green function in (M). In the equation of motion for this term we use the fact that 



n fi (S'c +S-c ) = l {Sln " ki ~ )c *t for (3 = -\- 

! " "'' ' ' ' c kl (Sy kl -(3S k a+) lor J 



and replace S k n k ^ — (3S k a k by its average in the Green function coming from Hq. 
Neglecting two-site correlation functions and commutators of 5nf, we obtain the equation 
of motion, for i ^ k, 



Ki< {S z k c k , + S k c kl ) ; c] T )) » (Sl4 x - (3S k at) £ f„ ((tfnfo ; c} t 

i 

"^'''V^cj,)), (8) 
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thus closing the system of equations for the scattering correction. 

We solve (|8]) for the Green function on the left-hand side, substitute the result into 
(0), and rearrange and sum over (3 to obtain, again for i ^ k, 



«<KH T ; c h)) t = ?r( £ )E f H «K<*t ; c 1t)) 6 > ( 9 ) 



where <?f (e) is the atomic-limit Green function presented in (f|). In the appendix of [p70f 
Hubbard solved this equation in terms of {(^nfjCq. ; )) e : 

Kl c *T : c iT» e = (E ((<KI c *T ; c] T )) e , (i ^ fc) (10) 
where Wj:f is defined by 

W£!(e) = g k J(e) - (Ha) 

^■( e ) = lv exp(ig -^ ) . (116) 

It is easy to check Hubbard's solution: from (|116|) it may be shown that g CT iife^ = 



5? — g'o-^y an d hence that g a Yl,i tiaWJf = WjZ — g a ^kj for i ^ k, and this can be used to 
verify that substituting ( fLOf) into (||) does indeed give a solution. 

As in Hubbard's case g l j is a zeroth order approximation to G l J in which both 
scattering and resonance broadening corrections are neglected, and following Hubbard 
we make the self-consistent replacement gl J \— > GV; this is the essential self-consistency 
step of the CPA. From fllPf ) the scattering correction term in (|5]) is given by 

E <« ((<K c *t ; c lT» e = ((k^t ; c lt» e ( 12 ) 

k 

where J a {e) = Y^ki^ikW^f (e)t^. A result that will be useful later is 

J CT (e) = e-S CT (e)-G (T (e)- 1 (13) 

which holds if the self-energy S CT (e) is local — as is expected for a CPA. This result may 
easily be established by making Fourier transforms and using G*J(e) = N" 1 X^exp(ifc ■ 
Rij)G^ a (e). We now also define the useful quantity 

E a (e) = e-J a (e). (14) 

It remains to find expressions for the unknown Green functions and Jf in (|5|); since 
these represent propagation as f- and j-spins respectively we neglect the commutators 
[nf,,H ] and [n"pif ] in the respective equations of motion: 

eSf a (e) « (Stnl) ^ + £t*G« (e)) + p jfc ((tf«5f)c fct ; c] T )) f - ^ a+ Tf Q (e) 

~ ((< ((^) 2 c iT + 5^r c .J ; C ] T )) £ (15) 
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(((<T 5 i c ^ + a ( s i a t) c ki) ; St 



j 



_1 («5f °( e ) - 5 Q _Tf °( e ) + 5(5 + l)Gf (, 

+^ ((< ((^) 2 c iT + 5^r Ci J ; c ] T )) e . (16) 

The second terms on the right-hand sides of these equations are the scattering 
corrections; the scattering correction in ( [15]) is of the same form as in ([|), but in (|16|) 
it is more complicated. It is not clear what the average 'zeroth order' Green function 
should be here, but we have chosen it to be —a(S^af)G k ^(e) as this makes the first 
terms on the right-hand sides of (|3p, (|1^), and (|IIJ) all of the same form; it turns out 
that this is necessary for the consistency of the approximation. The last Green function 
in ( |T5"D and ([16]) involves higher order spin operators and is in general unknown; this 
term must be treated using an approximation which is exact in the atomic limit. 

The scattering correction of (|T3D may be treated in the same way as that of (B); the 



only difference is that 5n^ is replaced by 5(n^S^), hence 

E^^K^jeJt)) ~ J T (6)((5«^)c iT ;c] T )) . (17) 



It is less obvious how to treat the scattering correction of flTq). Instead of giving 
a rigorous derivation we proceed by analogy and make the apparently reasonable 
approximation 

£*tk (((<T< S T c fc| + " (^ r ^ + ) c fcT ) ; cJ T )\ 
k 

« Me) ((n%Sr Cil ; + J T (c) ((a (S^f) c, T ; cJ T )) e (18a) 

= Jl(e)T? a (e) + a (S-af) Me)G^(e). (18b) 

It will be seen later that this is consistent with the rest of the derivation, leading for 
instance to a local self-energy as expected in a CPA. 

The system of equations of motion is now closed apart from the Green function 
((n"((S' J z ) 2 Cj-|. + S^S~c il );Cj-^)) e appearing in fll5|) and ([TB]). In general we will have to 
write more equations of motion to close the system. This is straightforward to do if 
we make approximations analogous to those used to obtain the equations of motion for 
5 T and T T . Then the Green functions («(5f) m c iT ; c] T )) e and {{nf(Sf ) m " 1 5"c ii ; cj T )) e 
and the expectations ((Sf) m n~,) and ((Sf ) m_1 S' l r cr l + ) for m = 1, . . . , 2S are brought 
into the system of equations. The resulting system of 4S + 1 equations per spin is 
complicated to solve self-consistently for large S however, so we will restrict ourselves 
to simple special cases in which we do not need any more equations of motion; in the 
next three subsections we will consider the paramagnetic state in zero magnetic field 
for arbitrary S, the case S = 1/2 for arbitrary magnetization, and the case of saturated 
ferromagnetism for arbitrary S. 

Our approximation scheme is now complete and is entirely self-consistent, i.e. all 
the expectations appearing in the system of equations can be obtained from the Green 
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functions in the system via the relation [13 

(BA) = - r -f(e-^)lm((A;B)) € (19) 

J-oo 7T 

where /(e) = 1/(1 + exp(/3e)) is the Fermi function and fi is the chemical potential. 
Also, since all the approximations made above have been in terms proportional to 
the approximation is exact in the atomic limit for all S and n as required. The system 
is assumed to be homogeneous, so we now drop the site indices of the expectations. 

3.1. Paramagnetism 

In the case of paramagnetism with zero field J-\ — Jj_. Hence after substituting fll7D and 
( P-8&D in ([TJI) and ([16]), respectively, we can add the resulting equations so that S^ a and 
T l i a occur only in the combination S^ a + T^ a : 

e (S ija (e) + T^(e)) = (S^ - aS'a^) (sy + ^> fc G^(e) - J(e) G^{e)) 

+ fa) - ^f) (S ija {e) + T ija {e)) - JS ^f l) G ija {e). (20) 

Since S^ a and T yQ also enter (|J) as S^ a + T^ a we have now closed the system. We 
substitute fll2"D and (p0|) into @ and rearrange to obtain 



(21) 



where G is defined by 



, v (E(e) + aJ/2)(n/2y + Ja/2(S-a) 
U Q ti( J E(e)-«J^/2)(E( e )+«J(S+l)/2)' 1 J 

E{e) = e — J(e), and (n/2) a = <5 Q _ + an/2. The spin symmetry of the paramagnetic 
state has been used to simplify the expectations appearing in (p2|). 



Taking Fourier transforms and solving for Gj; (^) becomes 

GAe) = i . (23) 

k J{e)+G{e)-^-t % 

Since in general G%(e) = (e — £g — S^(e)) -1 this implies that the self-energy is local, 
S?(e) = S(e), and so using expression fllTf ) for J(e) (|23|) becomes 

^r(e) = , (24) 

and for consistency in this equation we must have G = G. The existence of this simple 
self-consistency condition defining G is our main justification for the approximations 



made for the scattering correction of ([H]). If the approximations are changed we will 
not obtain a consistent expression defining a local self-energy. 
In the low density limit n — ► (|22|) reduces to 

G(e) = m ~ 3,2 (25) 

{ } (E(e) + JS/2)(E(e) - J(S + l)/2) ' 1 ] 



which is Kubo's [Q equation for the paramagnetic state Green function as required. 



3.2. The S = 1/2 case for arbitrary magnetization 

For S = 1/2 we have (S z ) 2 = 1/4 and SfS~[ = — 1/2 Sf for arbitrary magnetization, so 
the last Green function in QT5D and flrBj ) may be simplified: 

((< ((5f) 2 c jT + S*S- Cil ) ; c] T )) e = ~Gf*(e) - -if^e), (26) 
closing the system. Substituting the scattering corrections (|I2]), (0), and (|186fl into the 



equations of motion (||), (|15D, and (|16|) the system of equations reduces to 



E r G] ja = (npAf - - (S\ ja + If a ) (27a) 

7 To/ ■■ 

E^S\ 3a = (S'n^X] 3 - a - ^T T Jja (276) 

J Try ■■ 

E«T\ 3a = -a(S-a + )\] J - -G l f - ^Sf a , (27c) 



where 



A^'(e) = fiy +E**G«(e) _ J^G^e), (28) 

£/<j(e) = e — J«j(e), and E° = E a + Ja/4. These equations can be solved using a similar 
method to that of the previous subsection, with an analogous self-consistency condition 
occurring, and the local Green function is given by 



K) (E,(e)E?(e) - J 2 /8) - J/2 ((S^)E^(e) - a (S-^)E^(e) J 
h E^(e)(Ef(e)Ef(e)-jyA) 



The expectations (nf), (S z n"), and (S* a + ) may be calculated self-consistently using 
(WR) -<WR) and (|19D. 



If we set n = (|29|) reduces to 

r n _ ^(^)^(^)-^/8-J/2(g-)^ + (^) 

T( } " Ejw^mw^m • (30) 

For S = 1/2 the probabilities of a local spin being up or down are given by P(S Z = 
±1/2) = 1/2 ± (S z ), and this fact may be used to show that fl30|) is equal to Kubo's 
equation |14| for in the 5 = 1/2 case. Agreement with Kubo's approximation in 
the case of arbitrary magnetization is a much more stringent condition on our CPA 
than agreement just in the paramagnetic case, so this gives us confidence that our 
approximation is indeed a many-body extension of Kubo's. 

3.3. Saturated ferromagnetism for arbitrary S 

If we assume the existence of a saturated ferromagnetic state where all local spins 
and conduction electron spins are aligned parallel to the positive z-axis substantial 
simplification occurs. The expectations occurring in the |- and j-spin systems are given 
by 

n t f) = 6 a - (S z n a l ) = S5 a . (S~a + ) = (31a) 

nf) = 5 a - + an (-S z rf) = -S (5 a _ + an) /s + a~)=0, (31b) 
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and the Green functions G a , S a , and T a may be simplified using the relations 

G] j+ = T\ ]a = S\ ja = S5 a ^G { f (32 a) 

Tf + = Sf a = -SG l l a . (326) 

The higher order Green functions are given by 

< ({Slfc, + S*S~ Cll ) ; c] T )) = S 2 5 a ^ (33a) 



<(< ((5?)^ - SfSt^) ; c] A )) e = 5 2 Gf - S^_lf • (336) 

The system therefore closes without needing any more equations. After simplification 
the set of equations of motion (for both spin types) reduces to 

E ] + G} j = Af (34a) 

E i - ^) G* + = nAf (346) 

^-^)Gf" = (l-n)A?-^ (34c) 

(£ T + J(,S ~ ^ Tf = - J5 Gf (34d) 

with the other equations satisfied automatically. Here X l J and E a are defined as in the 
previous subsection. 

We can solve these equations in the same way as in the previous two subsections, 
with the usual self-consistency condition similar to (^4j) applying, and the local Green 
functions are given by 

G ' = E^TJST2 (35fl) 

r n | forn 

1 E i - JS/2 Ei- JS/2 - J 2 5/2 (£ T + J(S - l)/2)- 1 ' 1 ; 
The |-spin Green function here is just the free Green function shifted in energy, as 



might be expected. In the n — > limit these equations reduce to those of Kubo 14 



providing a further check on our approximation. It is shown in the next section that 
these Green functions are in fact not consistent with the initial assumption of a saturated 
ferromagnetic state. 



4. CPA spectral function 

In this section we study the DOS, first for the saturated ferromagnetic state and 
then for the paramagnetic state. We make the approximation of replacing the true 
DOS D c for a simple cubic lattice with the elliptic DOS, D e (e) = 2/(-nW 2 )^/W 2 - e 2 
where 2W is the bandwidth. For an elliptic DOS the free Green function is given by 
G (e) = 2/W 2 (e — \/e 2 — W 2 ). Since the self-energy E CT (e) is local the full Green function 
is given by G a (e) = G (e — E (T (e)), and it follows immediately that J a = W 2 /4G (T . 
The Green function equations therefore become algebraic in this approximation, greatly 
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simplifying the calculations. The elliptic DOS is also a fairly good approximation to the 
true DOS, unlike the Lorentzian DOS D\{e) = (W ' j V) / (W 2 + e 2 ) considered by Furukawa 
0|. In the infinite dimensional limit of a hypercubic lattice with nearest neighbour 
hopping — the scenario considered by Furukawa — it may be shown that the true DOS 
is a Gaussian if the hopping tij is scaled as the inverse square root of the number of 
dimensions [JT5| . With the scaling appropriate to three dimensions the Gaussian DOS is 
given by D g (e) = (3/ n) 1 / 2 exp[—3(e/W) 2 ] /W . These DOS's are compared in figure [TJ. 



1.0 




Energy (W) 



Figure 1. The DOS's D c , D e , D h and D g plotted in units of W. 



4-1. Saturated ferromagnetism 

We study the ferromagnetic state with complete spin alignment in the strong-coupling 
J — > oo limit which is most favourable to ferromagnetism, shifting the energy origin by 
— JS/2 in order to have the zero of energy in the lowest band. Equations ( |35a|) and 
( p56|) become 

G T = 1/£ T (36a) 

G ' = E^k- <"»> 

The DOS of this state is plotted in figure § for S = 1/2 and various n. Since the |- 
and |-spin DOS's are nonzero for the same range of energies it is clear that consistent 
saturated ferromagnetism does not occur within our approximation for any n. It may 
easily be shown that this is true for any finite S. Since strong ferromagnetism is expected 
to occur in the DE model for at least some parameter values this is a limitation of our 
approximation and suggests that our CPA may not be very good at low temperatures. 
Of course in common with the usual alloy CPA we cannot obtain a true Fermi liquid 
groundstate in our CPA — the imaginary part of the self-energy does not vanish at the 
Fermi surface at T = — so we only expect our CPA to describe the DE model well 
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at finite (but not necessarily large) temperatures where this is not a problem. The 
possibility of weak ferromagnetism occurring in our approximation will be discussed 
elsewhere; preliminary work on the magnetic susceptibility has already appeared in 



0.75 




Figure 2. The f-spin (solid curve, independent of n) and |-spin (dashed curves, 
n = 0, 0.25, 0.5, 0.75, DOS decreasing with increasing n) DOS's in the saturated 
ferromagnetic state for S = 1/2 and J = oo. Energy units of W are used. 



Note that the weights of the |- and j-spin bands here are 1 and (1 — n)/(2S + 1) 
respectively, so in the classical spin limit S — > oo there is no possibility of j-spin weight 
occurring near the Fermi level. In fact for S — > oo the saturated state is always stable 



against spin reversal but for finite S this is not always the case [}16 
4-2. Paramagnetism 

We now consider the zero field paramagnetic state. It may be shown that (S-a) — > nS/2 
as J — » oo, and (S ■ a) will be very near to this limit as long as JS <; 2W. We make this 
approximation in figure ^], in which the paramagnetic state DOS is plotted for 5 = 3/2 
and J = 4W for various n; this has the effect of removing the weak band centred on 
— J(S + l)/2. It may be seen that as n increases from the band near J(S + l)/2 is 
reduced in weight and a new band appears near JS/2, until at n = 1 no weight remains 
in the band near J(S+l)/2. The weight in the band near —JS/2 is (S+l—n/2)/ (2S+1) 
per spin, so if JS is sufficiently large to separate the bands (JS £ 2W) this band will 
just be filled at n = 1 producing a Mott insulator, as discussed in section |2|. 
We can understand figure |3| by expanding ( p2|) in partial fractions, 



G{e) 



1 v / (S + l)(n/2) a + (S-a) , S(n/2) a -(S-a) 
2S + l£l 



E(e) - aJS/2 



E(e) + aJ(S + l)/2 



(37) 



Comparing this equation with (W) in the paramagnetic state we see that G(e) = g(E(e)), 



13 




n=0.75 



0.50 



CO 
CO 

o 

>. 

*^ 

in 
c 
d> 
Q 



0.00 




n=1 



0.50 



re 
CO 
o 

</) 
c 


Q 



0.25 



0.00 




-0.5 1.5 

Energy 



Figure 3. The DOS in the paramagnetic state for S = 3/2, J = 4W, and 
n = 0, 0.25, 0.75 and 1. Energy units of W are used. 



and since E(e) — > e as — ^ it is clear how the bands centred on the peaks of the 
atomic limit spectral function arise as £y is switched on. 

In the strong-coupling limit J —>■ oo, which is taken with the energy origin shifted 
by -JS/2, (HD simplifies to 

G{e) ^ , (38) 

which corresponds to an elliptical band of weight (S + 1 — n/2)/ (2S + 1) and bandwidth 



2WJ(S + 1 — n/2)/(2S + 1). The band-narrowing factor is expected to depend upon 
the form of the bare DOS, as pointed out by Kubo |12|, but it is interesting to note 



that ours is the square root of the factor obtained by Brunton and Edwards ]16[ using 
a different method. 
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5. Comparison with the hole CPA 

In |H| Kubo and Ohata used a canonical transformation to derive an effective 
Hamiltonian H' for the J — > oo limit of the DE model. They then mapped H' onto an 
effective Hamiltonian H" for holes, in which conduction holes hop with spins aligned 
antiparallel to local spins (S + 1/2). In [T2j Kubo studied H" for n = 1, using a 



dynamical CPA similar to his n = CPA mentioned above. There is no obvious reason 
for one of these CPA's to be a better approximation than the other, so we would like 
our many-body CPA to agree with both of Kubo's approximations in the appropriate 
limits. In this section we set n = 1 and J = oo in our S = 1/2 Green function equation 



(f29|) and compare the result with Kubo's hole CPA. 

For S — 1/2 Kubo's n = 1 Green function is given by 

q _ 3/2 {Ei + gr) + {2E t + E^)(S Z + g) + 2(^ - £ T )(gV« ) 
T 2/3(2^ + ^(2^+^) 

where -E CT (e) = e — 3W /2 /8 G CT (e). If we set n = 1 and J = oo in fl2"9[) we obtain 

Gt = (gj + 3^ T )/4 + (gj + E,)(S* + a z )/2 + - E^)(S z a z ) _ 

E T (E T + Ei) 

For n = 1 and J = oo the relation (S 12 ) = (<j z ) holds, and we have used this to write 
both expectations as (S z + cr 2 )/2. Unfortunately the Green functions defined by (|3^) 
and ([K]) are not equal; this may be seen by considering the zero-field paramagnetic 
state where they reduce to 

respectively, and these equations correspond to bands with different widths. This 
appears to be a limitation of our CPA, which we would ideally like to interpolate 
continuously between Kubo's CPA's. A possible explanation for this disagreement is 
our neglect of resonance broadening corrections in section [5[ 

We could attempt to include resonance broadening corrections in our approximation 
by writing more equations of motion for these terms. It is however difficult to 
find approximations that give a closed set of equations for the new Green functions 
introduced without spoiling the self-consistency of our approximation. An alternative 
is to use an interpolation scheme containing arbitrary parameters which are chosen to 
yield as many correct moments of the spectral function as possible |18[ . 



We first note that an empirical modification of (0), of the form 



E^(e) = e- 



W_ 



2 



(42 a) 



3W 2 

E l (e) = e--£-G l (e), (42b) 

maps ( PD| ) for G-^ onto Kubo's equation This modification is of the same type as 

Hubbard's resonance broadening correction to the strong coupling (U — > oo) limit of his 
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model. This suggests that we should introduce an interpolation formula for E a which 



reduces to ( 42c ) and ( |42fcj) in the limits n = and n = 1, respectively. Unfortunately 



so far work along these lines has not proved successful. 
6. Resistivity formula 

We now derive a formula for the DC resistivity p of the paramagnetic state of our 
model, taking into account the cubic symmetry of the crystal and the local nature of 
our approximation. The Kubo formula [jn| states that for a small electric field uniform in 
space but oscillatory in time with frequency u>, E(r,t) = E exp(— iut), the conductivity 
tensor a^ u (uj) is given in terms of the current- current correlation function by 

/ \ r / / t t w /.ox 

where h — 1, Q is the unit cell volume, m and — e are the electron mass and charge 
respectively, and the retarded Green function is used. J is the electric current density 
operator defined for a homogeneous system by 

f =-m^ v ^ (44) 

ku 

where the velocity v% = Vgt^. In our case the conductivity is a real scalar so from 



a = —NQ lim Im 



(45) 



Now from (M 



e 2 



{{Jx\Jx)) w = 3(Nn) 2 E fy-0&{{ n & n P a >)) u > ( 46 ) 
^ ' kk'aa' 

so we need an approximation to the two-electron Green function ((n^ CT ; Tip^))^. 

Since the self-energy is independent of momentum a reasonable approximation is 
to assume that the irreducible vertex function is also independent of momentum. In 
infinite dimensions where the self-energy is rigorously local the momentum-dependent 
contribution of the irreducible vertex function vanishes. In this case the contribution 
of the vertex correction to a vanishes owing to the different parities of and t% in 
k, and we can evaluate the two-electron Green function ((n? ;nr in the bubble 
approximation, obtaining 

\j2 G kA iuJ rn + n)G nt7 {h} m + iu + n) (47a) 

P m 

= \l^l diE 7 \^ A kM {m) 

(3 J J V (iu m + p. - e) (iw m + iv + [i - 77) 
where the iu; m 's are odd Matsubara frequencies and iv is an even Matsubara 
frequency and the spectral representation of the one-electron Green function, A^ a (t) = 
— Im G^ a (e)/n, has been used. The sum over m can be evaluated using Cauchy's residue 
theorem, and following analytic continuation to the real axis and a shift in r\ we obtain 

»fc)) w = / de / d /^ (e) J^ + £) [f(e -ri-f(e + V - p)\ . (48) 
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In the paramagnetic state G is T- independent if we assume (S ■ a) = nS/2. The 
Fermi functions in ( fi"8|) give the conductivity a a weak T-dependence but we neglect this 
dependence and calculate at T = 0, considering our calculation to apply to the T > Tq 
state however. Equations ([45]), (fiBf), and (|48| ) then imply that 

2vre 2 



a 



3Ntt 

k 



( 49 ) 



Now for the simple cubic band = — 2t[cos(k x a) + cos(k y a) + cos(k z a)], where a is 
the lattice constant, and (^) can be simplified using Gauss' theorem: 

2T ' 2 E^fe) (50) 



a 



3Na 

k 



where = A%(ji) 2 . This is a legitimate definition since depends on k only 

through tj: in the local approximation. We introduce the cubic bare DOS D c (e), and 
(|5Cil) becomes 

lire 2 r 

o = — deeD c (e)<f>{e). (51) 



(52) 



3a 

Integrating A e (p) 2 with respect to e it may be shown that 

e + SY/i)-^ 1 , /e + S'(/^) - /i" 

tan 1 



2tt 2 



(E"(/i)) 2 + (e + S'(/i)-^) 2 S«(A*) V , 
where S'(yLt) and are the real and imaginary parts of £(//) respectively. 

If -D c (e) is replaced with the Gaussian D g (e) corresponding to an infinite dimensional 
approximation |I5[, (BTf) may be simplified by integrating by parts: 



2 TT/2 

a = / de J D g (e)^ e (//) 2 . (53) 



9a 

Note that this conductivity formula is of the same form as Furukawa's 



7. Resistivity calculations 



In this section we will use formulas fl5"TD , (|52j), and (1551 ) to calculate p = a -1 , making 
simple analytic approximations for D c (e). Note that in SI units ( pl| ) and (|53| ) must be 
divided by h. 



7.1. Elliptical DOS 

We first calculate p for the elliptical DOS D e (e). We take the strong-coupling limit 
J — » oo for simplicity; the precise value of J is unimportant in the strong-coupling 
regime JS > 2W 7 . p is plotted against n for various 5 in figure [|. It may be seen that 
p correctly diverges at n = and n = 1, but the size of p and dp/dn for n ~ 0.8 
are much too small to explain experiment, at least in some materials. For example 
Urushibara et al [IIJ, who measured p for La n Sr!_ n Mn03, found (for T well above Tq) 
that p ~ 20 mOcrn for n ~ 0.8 and that p drops by about an order of magnitude as n 
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is reduced from 0.85 to 0.7. We therefore agree with Millis et al [|[] that the bare DE 
model does not explain the paramagnetic state of the CMR materials. We also find that 
the T-dependence of p is always of the metallic dp/dT > form in our approximation, 
whereas Urushibara et al find a crossover between metallic behaviour at n — 0.7 to 
insulating behaviour at n — 0.85. It should be noted that p is independent of W here, 
and the size of J and S and the choice of resistivity formula flSlj ) or (|53|), using D e in 
both cases, does not have a large effect. 




Filling (n) 



Figure 4. The zero field paramagnetic state resistivity p (in mSlcm) versus filling n, 
plotted for J = oo, a = 5A, and S =1/2, 1, 3/2, 2, 5/2, and oo, p increasing with S. 
The elliptical DOS and resistivity formula fl5l| ) are used. 



7.#. Lorentzian DOS 

In order to make contact with Furukawa's work || we now consider the effects of 
approximating the true cubic DOS D c (e) with the Lorentzian D\(e), which is somewhat 
less realistic than the elliptical DOS approximation owing to its slowly decaying tails. 
Our equations are very simple in this case: it may be shown that J(e) = —iW and our 
equations for the Green function become explicit; switching on W here merely broadens 
the peaks of the atomic limit spectral function into Lorentzians with the same width 
parameter W as the bare DOS. Our approximation and Furukawa's become very similar 
in this case; in the classical S — > oo limit for example our local Green function is equal 
to Furukawa's. 

The slowly decaying form of the Lorentzian means that the precise value of J has 
more effect here than in the case of the elliptical DOS: the approximation (S-a) ~ nS/2 
only holds accurately for very large J and we only obtain a true insulator at half-filling 
for J = oo for instance. For simplicity however we take the limit J — > oo (with energy 
origin shifted by —JS/2) rather than use the finite value used by Furukawa. In the 
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paramagnetic state the local Green function G is then 

s + 1 _ n/2 i 

u 25+1 e + jW K J 

and the self-energy and chemical potential are 

E(e) = - f 0/0 2g t^ 1 ( e + iW ) ( 55fl ) 



/i = W tan 



2(5 + 1) -n 
7i f(4S + 3)n-2(S 



2 V 2(5+1) 




(556) 



If we use our finite-dimensional formula (|5"ID to calculate p with the Lorentzian DOS 
we must use a mixed approximation to the DOS, calculating using the Lorentzian 
approximation but using the elliptic approximation for D c (e), in order for the integral 
in Q5"Tp to converge. This problem does not arise when using the infinite dimensional 
conductivity formula ( p3| ) so D\ can be used throughout. The choice of conductivity 
formula flBT] ) (with mixed approximation) or ( |o*5| ) does not qualitatively alter the results, 
but since the calculation using flslj) has already been published |§ we use ( |53l) here. It 
also seems more satisfactory to use the same approximation for the DOS throughout. 
We believe conductivity formula ([51]) to be more realistic however, as it does not rely 
on the infinite dimensions approximation. 

p is plotted against n for various S in figure [5|. Note that the result is independent of 
W . As before we obtain insulating behaviour at n = and n — 1, but now p and dp/ On 
are of the same order of magnitude as experiment for n ~ 0.8, and our results are also 
consistent with those of Furukawa. Even better agreement can be obtained by reducing 
J, which has the effect of reducing p and dp/dn for n ~ 0.8. Furukawa's result — that 
the DE model's prediction for p is of the same order of magnitude as the experimental 
value — is therefore due to the unphysical approximation he used for the cubic DOS. 



8. Summary 

In addition to theoretical predictions of the properties of the DE model, experiments 
H also suggest that coupling to lattice degrees of freedom is a necessary ingredient in 
a model of CMR systems. However, we take the view that a good understanding of 
the bare DE model is necessary before more accurate and complicated models can be 
tackled. Accordingly in this paper we have presented a many-body CPA for the DE 
model and a compatible formula for the resistivity p; our main result is the reconciliation 
of the calculations of Millis et al Q and Furukawa of p in the paramagnetic state. 
We have confirmed that single-site scattering within the DE model is inadequate to 
describe the resistivity of CMR systems. 

Now the alloy CPA is known to give good results for a wide range of systems in 
which localization effects are unimportant. This should be the case here since, although 
some authors have claimed otherwise |2(| , localized states are expected to occur for the 
bare DE model only in the tails of the band; we therefore expect our formulae for G a to 
be good approximations. Nevertheless several problems exist with our approximation. 
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Figure 5. The zero field paramagnetic state resistivity p (in mf2cm) versus filling n, 
plotted for J = oo, a = 5A, and S =1/2, 1, 3/2, 2, 5/2, and oo, p increasing with S. 
The Lorentzian DOS and resistivity formula (|53"|) are used. 



Firstly, in common with the usual alloy CPA the imaginary part of the self-energy does 
not vanish at the Fermi level at T = so we do not obtain a true Fermi liquid; this should 
not be important in the high temperature paramagnetic regime that we have considered 
but the problem could perhaps be fixed following the approach of Edwards and Hertz 
to the Hubbard model f21fl . Secondly, we have shown that for finite S the saturated 
ferromagnetic state is never consistent in our approximation. Real CMR systems are 
believed to exhibit strong ferromagnetism in some regimes, and it is generally expected 
that the DE model has regions of strong ferromagnetism. Finally, although we recover 
Kubo's n = CPA in the low-density limit, we do not recover his strong-coupling n = 1 
CPA \ Y1\ at half-filling. This appears to be a shortcoming of our approximation, possibly 
due to our neglect of resonance broadening corrections. These problems and possible 
remedies will be discussed in a later paper in which we will concentrate on the magnetic 
susceptibility of our CPA. 
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